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The ionization energy based resistivity model with further confinements from spin-disorder scat- 
tering and polaronic effect is derived so as to capture the mechanism of both spin independent and 
spin-assisted charge transport in ferromagnets. The computed Tc-ossover below Tc and carrier den- 
sity in Gai_a;Mna:As (x = 6-7%) are 8-12 K and 10^^ cm~'^, identical with the experimental values 



of 10-12 K and 10^*- 
Lai-^iCa^^MnOs {x - 



10 cm respectively. The calculated charge densities for Mno.02Geo.9s and 
= 10-20%) are 10^^ cm~^ and lO^'^ cm"^ respectively. 
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1. Introduction 

Diluted magnetic semiconductors (DMS) have the 
tremendous potential for the development of spintronics 
and subsequently will lay the foundation to realize quan- 
tum computing. In order to exploit the spin assisted 
charge transport, one needs to understand the transport 
mechanism such as the variation of resistivity with tem- 
perature and doping in both paramagnetic and ferromag- 
netic phases. A wide variety of the magneto-electronic 
properties based on doping and Mn's valence state in 
manganites were reported to understand the transport 
mechanism(s) [l|,l2]- Among them, the influence of grain 
boundary as a barrier [^ , as a region of depleted Tc |j| 
and polaronic effect p on electrical properties were re- 
ported. Direct proportionality of H with Tc [g and 
resistivity with defects or substrate-film lattice incom- 
patibility [3, B IE U^ UM ^'^^ ^Iso regarded as equally 
important to determine the electrical properties of man- 
ganites. Furthermore, metallic conduction below Tc has 
been studied using double exchange mechanism (DEM) 
between s and d orbitals [13 and the displacement of 
hysteresis loop in field-cooled sample with an additional 
scenario of non-linear spin and charge fluctuations due 
to magnon 'l3|. Explanations in term of hopping elec- 
trons and DEM [lj| , and the influence of microstructural 
transition arises from ionic radius or valence state of Nd 
in Nda;Smi_j;Cao.8Mn03 [13 were also reported exten- 
sively. The effect of hydrostatic (external) pressure (P = 
0-^-15 kbar) and chemical doping (internal P) on metal- 
insulator transition of Pr-Ca,La-Sr-Mn03 have been re- 
ported as well. Ilq . 

Interestingly, Van Esch et al. [l3| have proposed multi- 
ple exchange interactions, which are ferromagnetic (FM) 
hole-hole and antiferromagnetic (AFM) Mn-hole interac- 
tions for DMS. These two effects, after neglecting the 
direct exchange between Mn-Mn (due to very diluted 
nature of DMS) are seem to be sufficient enough to de- 
scribe the temperature dependent magnetization curves 
{M{T)) accurately. However, even after inclusion of FM 



and AFM effects including the spin disorder scattering, 
the transport property in the FM phase is still not well 
understood. Unfortunately, this is also true for the case 
of metallic property below Tc in the well known and ex- 
tensively studied FM manganites as pointed out by Ma- 
hendiran et al. [l^. The resistivity {p{T)) above Tc for 
manganites is found to be in an activated form described 
by the equation |l8|. 



PiT > Tc) 
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(1) 



Ea is the activation energy, po and ks denote the resid- 
ual resistivity at T ^ Ea and Boltzmann constant re- 
spectively. In the FM phase, the influence oi M{T)/Mq is 
more pronounced than the electron-phonon (e-ph) contri- 
bution where the latter requires an overwhelmingly large 
coupling constant [l^]. Note that Mq is the magnitude 
of magnetization at K. Therefore, Mahendiran et al. 
have suggested that conventional mechanism namely, e- 
ph scattering has to be put aside so as to explain the p{T) 
for manganites below Tc- On the contrary, p{T) with e- 
ph involvement for DMS in the paramagnetic phase is 
given by [13 

,(r,T) C,+C,[e^p{Qn/T)-l]-' 

^ ^' fcsrin[l + exp((S„-^/)/fcsr)]' ^ ^ 

The term, C2/[exp(9_D/T) — l] takes care of the e- 
ph contribution. 0c , Ef, £",„, Ci and C2 represent 
the Debye temperature, Fermi level, mobility edge and 
numerical constants respectively. The p{T) in the FM 
phase based on the spin disorder scattering as derived by 
Tinbergen-Dekker is given by [l9j 



Psd{T < Tc) 



(m 



ejil 



Yo/2 



N{2Ei. 



Nl/2 



S{S + 1)-S^ 



TT{n,p)e^h'^ 
Mtd{T)'^ 



Mo 



S 



Mtd(T) 
Mo 



^ *""^ [ 2TS{S+l)Mo 



(3) 



N is the concentration of nearest neighbor ions (Mn's 
concentration) while {n,p) is the concentration of charge 
carriers (electrons or holes respectively) . m* ^ denotes ef- 
fective mass of electrons or holes, h = h/2-K, h = Planck 
constant, e is the charge of an electron, Ep and Jex are 
the Fermi and FM exchange interaction energies respec- 
tively while S is the spin quantum number. Equation (jS)) 
becomes equivalent to Kasuya 20] if one replaces the 
term, tanh[3TcAfTD(r)/2T5(5+l)Afo] with 1. Again, 
an accurate equation for the p{T) below Tc is still lack- 
ing since spin disorder scattering alone is insufficient as 
shown by Tinbergen and Dekker [l9j as well as reviewed 
by Ohno HJ. 

As a consequence, it is desirable to derive a formula 
that could describe the transport mechanism of ferromag- 
nets for the whole temperature range i.e., for both param- 
agnetic and FM phases and even at very low T. With this 
in mind, the Ej based Fermi-Dirac statistics (iFDS) and 
spin disorder scattering based resistivity models will be 
employed in order to derive p as a function of T, Ej and 
Mp (T, Mo ) . The consequences oi p{T,Ei, Mp{T, Mq ) ) 
that arises from the variation of T, Ej and Mp{T,Mo) 
are discussed in detail based on the experimental data 
reported by Van Esch et al. [ig, Mahendiran et al. |l8l | 
and Park et al. 22]. The Mn-cGei-^; FM semiconductor 
(FMS) is also accentuated here due to its promising prop- 
erties for device applications |22| where its gate voltage 
of ±0.5 V is compatible with the present Complemen- 
tary Metal-Oxide-Semiconductor (CMOS) and Ge's hole 
mobility (110.68 m^V^^s-i) is higher than GaAs (12.65 
m^V-^s-i) and Si (15.81 m^V-^s"^). Mn^Gei-^^'s re- 
sistivity is also semiconductor- like below Tq, which is 
more suitable than metallic Gai-^^Mn^jAs. Moreover, 
Mna;Gei_2, is also the simplest two-element system that 
can be utilized to evaluate the performance of the derived 
model consists of iFDS and Ma{T) (originates from tsd)- 
a = K (calculated from the Kasuya's spin disorder scat- 
tering model), p (calculated from the resistivity model), 
exp (determined experimentally). 



2. Resistivity model 

The total current in semiconducting ferromagnets with 
contributions from both paramagnetic and FM phases is 
J = X^y '^I'l ^ = 6''', se^ , h^, sh^ . For convenience, the 
spin-up, t denotes the direction of the magnetic field or a 
particular direction below Tc, while the spin-down, J, rep- 
resents any other directions. Note that the total energy 
(Kinetic -|- Magnetic), E\._^^ ^ ^k+m ^^^ *o energy 
level splitting below Tq- As such, the total current can 
be simplified as J = Jj- -|- Jj^ = Jg + Jse if the considered 
system is an n-type while J — Jh + Jsh if it is a p-type. 



Je and Jh are the spin independent charge current (elec- 
trons and holes respectively) in the paramagnetic phase 
whereas Jse and Jsh are the spin-assisted charge current 
in the FM phase. Thus the total resistivity {n or p-type) 
can be written as 



Pe.h + Pse.sh 



{n,p)e^Te 
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(4) 



Tsd represents the spin disorder scattering rate. The 
carrier density for the electrons and holes {n,p) based on 
iFDS are given by [Hil HI HE 113 
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iFDS is derived in a latter section in which, its applica- 
tions are well documented in the Refs. [HIm EIIm EtI 
UmillEl. Substituting I/t^ = AT"^ (due to electron- 
electron interaction), Eqs. Q and ©or © into Eq. Q), 
then one can arrive at 



.(T) 



ABexp[{Ej + EF)/kBT] 
AT^/^[Mp{T, Mo)]-i + BT-i/2 ■ 
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In which, A = [Ae,h/2e\mljy/^][27rhykBf/^ 
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[N{2EFf'\ml^f'^/TTh'']Jl,Mp{T,Mo). A.^n is the T 
independent electron-electron scattering rate constant. 
The El here takes care of the polaronic effect or more 
precisely, the electron-phonon interaction. The empirical 
function of the normalized magnetization is given by 



Mp(r,Afo) = l 



Mp{T) 
Mn ■ 



(8) 



Equation (jS)) is an empirical function that directly 
quantifies the influence of spin alignments in the FM 
phase on the transport properties of charge and spin 
carriers in accordance with Eq. (|2j). In other words. 



the only way to obtain 



Mp(T) 
A/n 



is through Eq. ©. In 



fact, Eq. ^ is used to calculate Mtd{T)/Mo and 



Mk{T)/Mo by writing S{S 

L2TS{S+l)Mo 



1) 



S(Mz^) tanh [ ^TSLslml ] = Mp{T,Mo) and 5(5 



Mo 



l)-^'(^^) -^(^^) = MpiT,Mo) respectively. 
Consequently, one can actually compare and analyze the 



Ma{T)/Mo (a = TD, K, p) calculated from Tinbergen- 
Dekker (TD), Kasuya (K) and Eq. 10 with the exper- 
imentally measured Mexp{T)/MQ. However, one has to 
switch to Eq. @ given below for the hole-doped strongly 
correlated paramagnetic semiconductors, which is again 
based on iFDS El El, 
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Ah is the T independent electron-electron scattering 
rate constant. Equation @ will be used to justify the 
importance of the term Jge even if the resistivity is 
semiconductor-like in the FM phase. Note that, m* = 
'^e ~ '^ft, ~ ("^e'^ft,)^ is used for convenience. If how- 
ever, m* ^ m)^, then one just has to use the relation, m* 
= m*m^/{m*^ -l-rn^J. Even in the usual consideration for 
the total conductivity, a — <J electron + crhoie, some alge- 
braic rearrangements can lead one to the relation, p{T) oc 
exp(£;//fcsr)/[exp(^i./fcBT) + exp(-^i./fci3r)], expos- 
ing the consistent effect of Ei on transport properties. 



3. Discussion 

3.1. Temperature- dependent resistivity curves 

Resistivity versus temperature measurement {p(T)) is 
the most simplest and effective method to study the 
transport properties. In free-electron metals, the p{T) 
curves are often exploited in order to deduce the T- 
dependence of the scattering rates namely, Te-e and 
Te-ph- Such behavior are well described by the Bloch- 
Grlineisen (BG) formula [31|, given by 



Pbg 



= Xt 



1287:m*kBT^ 



Bd/^T 



sinh^ 



-dx. (10) 



emission fine structure spectroscopy have exposed the po- 
laronic effect even at T < Tq or in the EMM phase |32| . 
Consequently, the charge density (n) in EMM phase is 
not T independent as one would anticipate for the free- 
electron metals (Fermi-gas). In addition, spin related 
mechanisms, like magnons and spin disorder scattering 
can be correctly represented with the normalized mag- 
netization function, M[T,Mq). It is quite common to 
employ Matthiessen's rule (r^^ = "^iT^^) as opposed 
to the total current rule (r = ^^ r^) and write the resis- 
tivity below Tc in the form of 



P{T) 



Po 



Y.'^^T" 



(11) 



The i here indicates the types of T-dependent scatter- 
ing rate that contribute to the resistivity and A is a T 
independent constant, po is the T- independent scattering 
rate that originates from the impurities as T — > K. The 
critical issue here is not about the Matthiessen's rule, but 
on the validity of Eq. (|ll|l in non free-electronic phase. 
Importantly, the T-dependent structure of Eq. (|l()|l is 
equivalent to Eq. (|ll|l that actually have enabled one 
to reliably calculate t~X{T) and T;:]:py^{T) as Ae-eT"^ 
and Ae-phT^^^ respectively. Equation ((TT1) is extremely 
popular and it is applied indiscriminately to determine 
the T-dependence of a wide variety of scattering rates 
in EMM phase, while the correctness of such determi- 
nation is still unclear and varies from one researcher 
to another \\A Is.li yj| . It is important to realize that 
only a free-electronic EMM phase at T < Tc will jus- 
tify the analysis based on Eq. Hll() . The influence of 
polaronic effect and magnetization function (the varia- 
tion of M{T)/Mo with T) reinforces the T-dependence 
of charge density, which point towards the inapplicabil- 
ity of Eq. (HTJ in EMM phase. 



3.2. Gai-xMn^As 



Xtr = electron-phonon coupling constant, n = free elec- 
trons concentration. The approximation of Te-e(T) and 
Te-ph{T) using Eq. (|10|) is valid basically because there 
are no other parameters that vary with T, apart from 
the said scattering rates. In fact, by utilizing the BG 
formula, one can reliably estimate that Te-e{T) oc T"^ 
while Te^ph{T) oc y-3^-5 for any experimentally viable 

On the other hand, the metallic phenomenon observed 
in the ferromagnetic-metallic (EMM) phase below Tc 
in ferromagnets (FM) cannot be characterized as Fermi 
gas. Therefore, it is rather incorrect to extract Te-e(T), 
Te^ph{T) and Traagnons{T) from thc p{T < Tc) curves 
in FM. Experimental evidences based on the photoe- 
mission. X-ray emission/absorption and extended X-ray 



The resistivity measurements |l7j and its fittings based 
on Eqs. and © are shown in Eig.Qa) and b) respec- 
tively for Gai_a;Mna;As. One needs two fitting param- 
eters {A and Ej) for p{T > Tc) and another two {B 
and Mp{T,Mo)) for p{T < Tc). All the fitting parame- 
ters are listed in Table HI Note that S' = 1 and 5/2 are 
employed for the fittings of Mk(T)/Mq while Tc and 
Tcrossover = Ter were determined from the experimen- 
tal resistivity curves. The deviation of Mk{T)/Mq from 
the Mexp/Mo increases with 5" from 1 -^ 5/2. The p{T) 
is found to increase with x from 0.060 to 0.070 due to 
the mechanism proposed by Van Esch et al. [ijl, |23 and 
Ando et al. Sfij. They proposed that neutral Mn'^+ ac- 
ceptors that contribute to magnetic properties could be 
compensated by As, where for a higher concentration of 



Mn, instead of replacing Ga it will form a six-fold coordi- 
nated centers with As (Mn^^'') [Hill Hi] • These centers 
will eventually reduce the magnitude of ferromagnetism 
(FM) in DMS due to the loss of spin-spin interaction be- 
tween Mn(3d^) and h. iFDS based resistivity models also 
predicts that if one assumes Mn^+ {Ej = 1113 kJmol"^) 
or Mn3+ {Ei = 1825 kJmol-i) substitutes Ga3+ {Ej = 
1840 kJmol"-^), then p{T) should further decrease with 
X, which is not the case here. Thus, iFDS also suggests 
that Mn^+ or Mn^+ do not substitute Ga^+. Interest- 
ingly, the Tc,s observed in Gao.94oMno.o6oAs (annealed: 
370°C) and Gao.gaoMno.oyoAs (as grown) are 10 K and 
12 K, which are identical with the calculated values of 8 
K and 12 K respectively. Note here that Ej + Ep = Tcr- 
The calculated carrier density using Ei + Ep (8, 12 K), 
m*j^ = rest mass and Eq. ^ is 2.4 x 10^^ cm^"^. Below 
Tc, spin alignments enhance the contribution from Jse 
and reduces the exponential increase of p{T). This re- 
duction in p{T) is as a result of dominating J^g and the 
small magnitude oiEi -\- Ep (8 K,12 K), consequently its 
effect only comes at low T as clearly shown in Fig. ^a). 
The Gao.g3oMno.o7oAs samples after annealing at 370 °C 
and 390 °C do not indicate any FM [l^ (Fig.db)). Thus 
the fittings are carried out with Eq. Q that only requires 
two fitting parameters namely, A and Ej + Ep since J^e 
= (there is no observable Tc) and/or dMaiT)/ModT 
— {Mp{T,Mo) — constant). The exponential increase 
in Fig.mb) for p{T) is due to Ej + Ep from Eq. ® with 
zilch Jse contribution. 

Figure ^c) and d) indicate the normalized magnetiza- 
tion, Ma{T)/Mo. Note that Mp,td.k{T)/Mq is a fitting 
parameter that has been varied accordingly to fit p{T < 
Tc). Mp{T,Mo) is used to calculate Mp,td.k{T)/Mq 
with S = \. M p,TD .k{T) I Mq is also compared with 
the experimentally determined [IJl Mexp{T)/AIo as de- 
picted in Fig. n d). One can easily notice the in- 
equality, Mtd{T)/Mo > Mk{T)/Mo > Mp{T)/Mo > 
Mexp{T)/Mo from Fig.Hc) and d). As such, Mp{T)/Mo 
from Eq. (O is the best fit for the experimentally mea- 
sured Mexp{T)/Mo. However, Mp{T)/Mo is still larger 
than Mexp(T)/Mo, because resistivity measures only the 
path with relatively lowest Ej and with easily aligned 
spins that complies with the principle of least action. On 
the contrary, the magnetization measurement quantifies 
the average of all the spins' alignments. 



is approximately 10^^ cm ^. In the presence of H = 6 
T, Ei + Ep is computed as 0.0776 eV for x = 0.2 that 
subsequently leads to p = 10^* cm~^. It is proposed that 
the activated behavior for p{T > Tc) is due to electron- 
phonon interaction or rather due to the polaronic effect 
{Ei) |23|. The fittings are shown in Fig. [21a) and b) while 
its fitting parameters are listed in Table Q] Theoreti- 
cally [13, Ca2+ {Ei = 868 kJuior^) < La3+ {Ei = 1152 
kJmol"""^), therefore p{T) is expected to decrease with 
Ca^"*" doping significantly. On the contrary, only a small 
difference of Ei + Ep between a; = 0.1 (0.12 eV) and 
0.2 (0.11 eV) is observed due to Mn''+'s compensation 
effect. The quantity, Mn^+ increased 6% from x = 0.1 
(19%) to 0.2 (25%) O. Ideally, the difference of Ej be- 
tween Ca2+ and La3+ is 1152 - 868 = 284 kJmor^. As 
a result of compensation with 6% Mn^+^^-t- ^^^ ^ 4940 
kJmol-i), the actual difference is only 284 - [0.81(1825) 
+ 0.19(4940) - 0.75(1825) - 0.25(4940)] = 97 kJmor^. 
This calculation simply exposes the compensation effect. 
At 6 T, Lao.8Cao.2Mn03 indicates a much lower resis- 
tivity (Fig. 121 b)). The result that larger H giving rise 
to overall conductivity is due to relatively large amount 
of aligned spins that eventually gives rise to Jse- Hence, 
Ei + Ep at 6 T {78 meV) is less than Ei + Ep a,t 
T (112 meV). Figure (21 c) and d) depict the calculated 
M„(T)/Mo with 5 = 1 and Mexp{T)/Mo for x = 0.2 re- 
spectively. The calculated Mtd{T)/Mo is dropped for 
Lai_2;Caa;Mn03 since Mk{T)/Mq seems to be a better 
approximation than Mtd{T)/Mq as indicated in Fig. Q] 
c) and d). The discrepancy between Mp{T)/Mo and 
Mexp{T) / Mq still exists even though Eq. (0) reproduces 
p{T) at all T range accurately. Again, this incompat- 
ibility is due to the principle of least action as stated 
earlier. In addition, the manganites' charge transport 
mechanism below Tc is also in accordance with Eq. (O 
because the term, Mp{T, Mq) handles the exchange in- 
teractions' complexities separately for DMS and man- 
ganites. For example, one can clearly notice the different 
type of discrepancies between DMS and manganites by 
comparing the empirical function, Ma{T)/Mo {a = p, 
exp) between Fig.Qd) and Fig. (21 d). Hence, Eq. (01 is 
suitable for both types of ferromagnets, be it diluted or 
concentrated. 



3.4. Mn^Gei- 



3.3. Lai-xCoxMnOs 

Mahendiran et al. [13 discussed p{T < Tc) with re- 
spect to Eq. (0) and obtained the activation energy, Ea 
= 0.16 eV for x = 0.1 and 0.2 of Lai_a;Caa;Mn03 sam- 
ples at T. Using Eq. however, Ei + Ep for the 
former and latter samples are calculated to be 0.12 and 
0.11 eV respectively. The calculated carrier density using 
Ei + Ep (0.12, 0.11 eV), m^ = rest mass and Eq. ® 



The Mn2;Gei_a; FMS with homogeneous Mn concen- 
tration has been grown using low-T MBE technique [23| . 
The Mna;Gei_a; was found to be a p-type with carrier den- 
sity in the order of lO^^- 10^° cm'^ for 0.006 <x< 0.035 
as measured by Park et al. |22| • Both the resistivity mea- 
surements I22I and its fittings based on Eq. Q are shown 
in Fig.Oa). Here, Ei + Ep, A and B have been floated 
while Mp{T,Mo) is constrained to reduce with T in or- 
der to fit the experimental p{T, x = 0.02). The absence of 



the Curie- Weiss law in the p(T, x — 0.02) curve is due to 
insuffi-cient number of aligned spins that eventually leads 
to a relatively small J^e, which in turn, is not able to pro- 
duce the metallic conduction below Tq. This scenario is 
also in accordance with the measured Mexp(T) where, 
only 1.4-1.9 /is/Mn atom contributes to ferromagnetism 
as compared with the ideal value of 3.0 /is/Mn atom. In 
other words, only 45-60% of Mn ions are magnetically ac- 
tive 113. It is found that Ei + Ep = 15K from the p{T) 
fitting for Mno.02Geo.98. Subsequently, one will be able 
to calculate the hole concentration as 2.38 x 10^^ cm~^ 
using Eq. |(BJ and m'^ = rest mass, which is remarkably 
in the vicinity of the experimental value |23|, 10^^ - 10^° 
cm~'^. Interestingly, the semiconductor- like behavior of 
p{T, X = 0.02) below Tq is not exponentially driven as 
the value of Ei + Ep is very small (15 K) to be able to 
contribute significantly above 15 K. Rather, it is the T- 
dependence of Eq. lO determines p{T, x — 0.02) below 
Tc- To see this clearly, one can actually approximate the 
experimental p(r, x — 0.02) with a mathematical func- 
tion given by p = -21.711 x InT + 148.47 (not shown). 
In this computation, the InT behavior is the approxi- 
mate version for the T-dependence of Eq. Q. Another 
obvious proof is the inability of Eq. © to represent the 
experimental p{T,x — 0.02). The plot using Eq. is 
also shown in Fig. |31a) with A^ ~ 1.8 and Ej + Ep ~ 
80 K that eventually give p ~ 1.92 x 10^^ cni^-^. How- 
ever, in the absence of Jse term, Eq. J^ is inadequate to 
capture the T-dependence of p{T,x = 0.02) in the FM 
phase. 

The pronounced effect of Mp{T, Mq) can be noticed 
by comparing the calculated plots between Eq. ijjj) and 
Eq. lO with additional constraint, dMp{T)/dT = as 
indicated in Fig. Oa). Recall that Mp(T,Mo) is varied 
with T to fit the experimental p(T, x = 0.02) in compli- 
ance with Eq. ^. Furthermore, p(T) is found |22| to 
decrease with x from 0.016 to 0.02 while p(T,x = 0.02) 
remains identical with p(T, x = 0.033). This type of tran- 
sition can be readily evaluated with Eq. O. Firstly, 
notice the large increase in room temperature p from 
10^^ cm~^ (upper limit) for pure Ge to 10^^ cm~^ (lower 
limit) for a mere 2% Mn substituted Mno.o2Geo98, which 
gives rise to a rapid decrease of p{T,x). The average 
Eis for Mn2+, Mn3+ and Ge'*+ are computed as 1113, 
1825 and 2503 kJmol"^ respectively. According to iFDS, 
Mn substitution into Ge sites will reduce the magnitude 
of p{T) since Ej (Ge'*+) > Ej (Mn3+) > Ei (Mn2+), 
regardless of dM{T)/dT = or dM(T)/dT £0. Such 
behavior has been observed experimentally [22| where, 
p(T,x = 0.009) > p{T,x = 0.016) > p{T,x = 0.02). 
This scenario indicates that the holes from Mn^+'^+ is 
kinetically favorable than the intrinsic holes from Ge^+. 
It is also found experimentally that Prt{x = 0.009) < 
Prt{x — 0.016) as a result of the variation of T indepen- 
dent scattering rate constants {A and B). Surprisingly 
however, p{T,x = 0.02) ~ p{T,x = 0.033), which sug- 



gests that A, B and T-dependence of Mexp{T) / Mq are 
identical. Using iFDS, one should get p{T,x = 0.02) > 
p{T,x^ 0.033) and [A,S],=o.02 > [A, S],=o.033. Mean- 
ing, the additional Mn substitution (0.033 — 0.02 = 
0.013) may not have substituted Ge, instead it could 
have formed a well segregated impurity phase that even- 
tually contributes to the higher magnitudes of A and B 
{[A, B]x=o.o2 « [A, B]x=o.o33), and consequently does not 
interfere with the Mexp{T)/MQ. Notice that the forma- 
tion of impurity phase is quite common in any system, 
including Gai-ajMn^^As DMS with strictly limited Mn 
solubility. On the other hand, the normalized magneti- 
zation, MK.p,expiT)/Mo for Mno.02Geo.98 have been plot- 
ted in Fig.Olb). One can notice the relation, Mk{T)/Mq 
> Mp{T)/Mo > M,xp{T)/Mo from Fig. Elb). Again, 
Mp(T)/Mo > MexpiT)/Mo is due to the abihty of both 
Je and Jse to follow the easiest path. Additionally, 
the T-dependence oi Mexp{T) is similar to Gai-ajMn^^As 
rather than the well established manganite ferromagnets, 
which reveals the possibility of multiple exchange inter- 
actions [13, U2, 1^ ■ All the values of Ei discussed above 
were averaged in accordance with Ei[X^'^] = X]^=i ■^• 
Prior to averaging, the 1**, 2"^, 3'"'* and 4*'* ionization 
energies for all the elements mentioned above were taken 
fromRef. [13 • 

4. iFDS 

Here, the Lagrange multipliers, degeneracy and the to- 
tal energy requirement associated with Ej in iFDS is dis- 
cussed in detail. Both FDS and iFDS arc for the half- 
integral spin particles such as electrons and holes. Its to- 
tal wave function, '5 has to be antisymmetric in order to 
satisfy quantum-mechanical symmetry requirement. Un- 
der such condition, interchange of any 2 particles {A and 
B) of different states, ipi and V'j (i 7^ i) will result in 
change of sign, hence the wave function for Fermions is 
in the form of 



*.,j(Ca,Cb) = MCa)^j{Cb) - V^,(Cb)V,(Ca). (12) 

The negative sign in Eq. (|12|) that fulfils antisymmet- 
ric requirement is actually due to one of the eigenvalue 
of exchange operator [38J, P = — 1. The other eigen- 
value, P = -1-1 is for Bosons. Ca and Cb denote all the 
necessary cartesian coordinates of the particles A and B 
respectively. Equation ((T^ is nothing but Pauli's exclu- 
sion principle. The one-particle energies Ei, E2, ..., Em 
for the corresponding one-particle quantum states gi, 92, 
..., qm can be rewritten as {Eis ± Ej)i, {Eis ± £■7)2, •■•, 
{Eis ± Ei)m- Note here that E.^s = Einniai state- It is 
also important to realize that Eis + Ej = Eeiectrons and 
Eis - Ej = Ehoies- Subsequently, the latter {E,s ± -E/)* 
version where i = 1, 2, ..., m with Ei as an additional 
inclusion will be used to derive iFDS and its Lagrange 



multipliers. This ±£'7 is inserted carefully to justify that 
an electron to occupy a higher state N from initial state 
M is more probable than from initial state L if condition 
Ej{M) < Ej{L) at certain T is satisfied. As for a hole 
to occupy a lower state M from initial state N is more 
probable than to occupy state L if the same condition 
above is satisfied. Eis is the energy of a particle in a 
given system at a certain initial state and ranges from 
+CX) to for electrons and to — cx) for holes. In con- 
trast, standard FDS only requires Ei {i — 1, 2, ..., m) as 
the energy of a particle at a certain state. 

Denoting n as the total number of particles with ni 
particles with energy (Eis ± Ej)i, 712 particles with en- 
ergy (Eis ± Ej)2 and so on implies that n — ni + n2 
+ ... + n,„. As a consequence, the number of ways for 
gi quantum states to be arranged among ni particles is 
given as 



dF 

- — =\n{qi - Tii) -In n.i + fi + X{Eis ±Ei)i 

an., 

= 0. 



(19) 



Thus, the Fermi-Dirac statistics based on ionization 
energy is simply given by 



(20) 



qi eyip[^i + X{Eis±Ei),] + l' 



Importantly, the total energy, E in iFDS can be ob- 
tained from Eq. 1)16(1 . which is 



E ^y(E,s±Ei),n, 



P{ni,qi) 



91 ! 



ni\{qi -ni)!' 



(13) 



Now it is easy to enumerate the total number of ways 
for q quantum states [q = qi + q2 + ■■■ + f?™) to be 
arranged among n particles, which is 



P{n,q)^\l 



^ nj!(gj -ni)\' 



(14) 



The most probable configuration at certain T can be 
obtained by maximizing P{n^ q) subject to the restrictive 
conditions 



V -^ ru2 

-^ 2m 



[k^±k?], 



fi2 fi2 

= — [k? ik^l =— k2. (21) 

k/ = Konized state, and the ± sign is solely to indi- 
cate that the energy corresponds to electrons is ^ +00 
while -^ —00 is for the holes, which satisfy the particle- 
hole symmetry. Consequently, Eq. (|21l) also implies that 
iFDS does not violate the degeneracy requirements. By 
utiHzing Eq. ^^ and taking exp[^ + \{E ± Ej)] > 1, 
one can arrive at the probability function for electrons in 
an explicit form as 



y Ui — n, y^ drii = 0. 



(15) 



/e(kis) = exp 



— /i — A 



2u2 



2m 



Ei 



(22) 



^(£,, ± Ei\iM = E, ^(£». ± Ei),(nM = 0.(16) 

i i 

The method of Lagrange multipliers [S^] can be em- 
ployed to maximize Eq. I|14|l . Hence, a new function, 
F{xi,X2, ...^i,\, ...) = / + m/i + •^/2 +••• is introduced 
and all its derivatives are set to zero 



dxn ' d^ ' 9A 



(17) 



As such, one can let the new function in the form of 



F = lnP + nJ2drM + A^(£,, ± Ei),dn,. (18) 



Similarly, the probability function for the holes is given 



by 



fhO^is) = exp 



fi + X 



2m 



~Ei 



(23) 



The parameters fi and A arc the Lagrange multipliers. 
h = /i/27r, h — Planck constant and m is the charge 
carriers' mass. Note that E has been substituted with 
h^'k^/2m. In the standard FDS, Eqs. ^ and ^ are 
simply given by, /e(k) = eyip[-^i-\{h^'k^ I2m)\ and A(k) 
= exp[/x -|- A(fi^k^/2m)]. Equation ifT^ can be rewritten 
by employing the 3D density of states' (DOS) derivative, 
dn = l^k|^dkis/27r^ Eqs. ^^ and ^, that eventually 
give 



After applying Stirling's approximation, dF/drii can 
be written as 



V 
2^' 



exp 



' 2m 



dk 



V 



re ^ J K exp 





V 
2^' 



= ^e"--^^^ /k^exp 



2m 2m 



dk. 



A- 



u 

P= 2^^'^-^'VkLexp 



2m 



21,2 



n^k 



2?Ti 



dk,,, (24) 



dk,,. (25) 



The respective solutions for Eqs. p4(l and (|25|l are 



fJ, + XEi 



In 



/i — XEj = In 



n 


f2TTxn^y^' 


V" 


\ m J 


P ( 


'2nXh^V^'] 


V \ 


^ m ) ^ 



(26) 



(27) 



Note that Eqs. (|2S1) and |53) simply imply that 
IJ.e{iFDS) = n{T = 0) + A^;/ and HhiiFDS) = /x(T = 0) 
— XEi. In fact, ii{FDS) need to be varied accordingly 
with doping, on the other hand, iFDS captures the same 
variation due to doping with XEj in which, /z(T = 0) is 
fixed to be a constant (independent of T and doping). 
Furthermore, using Eq. p6|) . one can obtain 



E 



AmTr"^ 

VfP_ 
Armr'^ 

Amn'^ 

3V 

2A' 

3V 



oo 

-KFDS) /j^^exp 



oo 

--(^=<^)|ktcxp 



oo 
g-M(T=0)-ABj f^ 

fj.{T=0)-\Ei 



ft2k2 



2m 



dk 



_A^^-A^ 
2m 2r?T, 

h^kl 



dki. 



:exp 



3/2 



A- 



2m 



dk,, 



2A 



-A'(-F-DS) 



27rAfi2 

-|3/2 



27rAft2 



(28) 
(29) 



Again, Eq. H28|l being equal to Eq. H29|l enable one 
to surmise that the total energy considered in FDS and 
iFDS is exactly the same. Quantitative comparison be- 
tween Eq. (|28|) and with the energy of a 3D ideal gas, 
E = 3nkBT/2, after substituting Eq. ^ into Eq. ^ 
will enable one to determine A. It is found that A remains 
the same as l/ksT. Basically, at constant temperature 
(T > 0), FDS predicts the distribution spectrum if E is 
varied, relying on external inputs such as band gap (Eg) 
and/or Fermi level (EpiT)). On the other hand, iFDS 



needs only Ej as an external input to predict the vari- 
ation of E, without relying on Eg and/or Ef{T) at all, 
and subsequently its distribution spectrum can be ob- 
tained as well. Notice that Ep comes into iFDS as Ep 
= constant, independent of T and doping. Ep denotes 
the Fermi level at K. In fact, iFDS and FDS take dif- 
ferent approach in term of energy levels and Fermions' 
excitations to arrive at the same distribution spectrum. 
In simple words, iFDS is new in a sense that it gives one 
an alternative route to obtain the Fermions' distribution 
spectrum in which, FDS needs Ep{T) and/or Eg while 
iFDS needs only Ej to arrive at the same distribution 
spectrum. Hence, based on the accuracy of these input 
parameters, one can choose either FDS or iFDS to be 
used for one's theoretical models. 



5. Conclusions 

In conclusion, the transport properties of 
Gai_a;Mn2;As, Lai-ajCaxMnOa and Mua^Gci-a; can 
be characterized with a model consists of ionization 
energy based Fermi-Dirac statistics coupled with spin 
disorder scattering mechanism. This model has been 
able to explain the evolution of resistivity's curves with 
respect to temperature and Mn doping. The arguments 
for the incompatibility between the calculated and 
experimentally determined normalized magnetization 
is based on the total current's tendency to obey the 
principle of least action. The validity of Ej + Ep 
and Mp{T, Mq) have been justified quantitatively by 
computing p and Mp(T)/MQ respectively, which are 
in excellent agreement with the experimental results. 
However, the magnitudes of A and B are not diag- 
nosed due to unknown reliable values of Ah, J ex and 
Ep. To this end, the variation of hole mobilities and 
dielectric constant with doping, the influence of multiple 
exchange interaction and energy gap above Tc should 
be investigated experimentally. 
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